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' Abstract. The rare earth pyrochlore magnet Yb2Ti207 is among a handful of 

. materials that apparently exhibit no long range order down to the lowest explored 

temperatures and well below the Curie- Weiss temperature. Paramagnetic neutron 
scattering on a single crystal sample has revealed the presence of anisotropic 
correlations and recent work has led to the proposal of a detailed microscopic 
' Hamiltonian for this material involving significantly anisotropic exchange. In this 

£h ■ article, we compute the local sublattice susceptibility of Yb2Ti2 07 from the proposed 

model and compare with the measurements of Cao and coworkers (HJ, finding quite 
good agreement. In contrast, a model with only isotropic exchange and long range 
magnetostatic dipoles gives rise to a local susceptiblity that is inconsistent with the 
. data. 
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• . 1. Introduction 

CN 
O 

The pyrochlore lattice of corner-sharing tetrahedra is a classic example of a lattice prone 
to geometric frustration, the inability of interacting magnetic moments on such a lattice 
to simultaneously minimize all of the pairwise interactions in the system. In the rare 
earth oxides, A2B2O7, where A is a trivalent rare earth ion (Ho, Dy, Tb, Gd, Yb) or 
yttrium (Y) and B is a tetravalent transition metal ion (Ti, Sn, Mo, Mn), both A and 
B reside on distinct interpenetrating pyrochlore lattices. Such systems have been the 
subject of intense research for the past 10 years as they are known to manifest a wide 
variety of exciting collective phenomena [T] including spin ices (Dy2Ti20"7, Ho2Ti207) 
[2], spin liquids (Tb 2 Ti 2 7 ) [3], spin glasses (Y 2 Mo 2 7 ) jl], and LRO with persistent 
spin dynamics (Gd2Sn 2 07, Er 2 Ti 2 07) (SlIEj. 

To better understand the many complex phenomena observed in magnetic 
pyrochlore oxides, knowledge of the microscopic interactions present in these materials 
is necessary. The pyrochlore lattice can be described as an FCC Bravais lattice with 
a tetrahedral basis at each FCC site. The magnetic rare earth {A 3+ ) ions reside 
on the pyrochlore lattice and are subject to a crystal field (CF) which introduces a 
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magnetic anisotropy with respect to the local (111) axes that is reflected in the magnetic 
susceptibility on each of the four sublattices. Recently, determination of the local 
susceptibility, Xa, via polarized neutron scattering has been brought to bear on the rare 
earth pyrochlores including Tb 2 Ti 2 7 , Ho 2 Ti 2 7 , Er 2 Ti 2 7 , and Yb 2 Ti 2 7 [8]. These 
measurements have been used to analyse the single ion anisotropy and the nature of the 
magnetic interactions [3 [9] . In particular, despite concerns about some details of the 
method with which the theoretical analysis of Xa is carried out in Refs. [H |9] , which we 
discuss later on in the Appendix, both of these works arrive at the tentative conclusion 
that in most of these materials, isotropic exchange interactions are inadequate to account 
for the Xa measurements and that, instead, the interactions appear to some extent 
anisotropic. The presence of anisotropic exchange in the pyrochlore oxides provides 
a further clue to the vehicle for the richness of phenomena observed in rare earth 
frustrated magnets p. Indeed, recent theoretical work has exposed large anisotropic 
superexchange in pyrochlore oxide materials [TU| [TTj. 

In this paper, we focus on local susceptibility measurements on one material, 
Yb 2 Ti 2 7 [8], which is remarkable for exhibiting a phase transition at T c w 240 mK, 
albeit with no signs of long range magnetic order at lower temperatures [12]. The 
nature of the low temperature phase in Yb 2 Ti 2 7 is currently somewhat controversial. 
Reference [13] reports the onset of hysteresis in the vicinity of 240 mK and the presence 
of magnetic Bragg peaks at lower temperatures. Other studies, however, have found 
no evidence of magnetic reflections [121 O EH] i n zero applied field. Reference [12] 
presents muon spin relaxation data revealing that the spin fluctutation rate drops 
sharply at about 240 mK, but remains nonzero and temperature independent well below 
T c , perhaps indicating a spin gas to spin liquid type phase transition. Peculiarly, the 
collective paramagnetic state exhibits strongly anisotropic rod-like correlations in the 
diffuse neutron scattering which have been conjectured to arise from a weak decoupling 
of the kagome layers in the pyrochlore lattice [15]. We have investigated the origin 
of the correlations seen in neutron scattering by using a random phase approximation 
calculation of the paramagnetic neutron scattering intensity pattern and found a set of 
interactions that produce a theoretical scattering pattern that is in very good agreement 
with the experimental data [16]. The foremost characteristic of these interactions is 
that they are strongly anisotropic. The purposes of the present work are (i) to provide 
a comparison of the predictions of our microscopic Hamiltonian [16] with the local 
susceptibility (xa) data of Ref. [8] and (ii) to provide a discussion of the sensitivity of Xa 
to the nature of the interactions and compare with the results of Refs. [8] and [9]. We 
show that our model can successfully account for the observed temperature dependence 
of Xa without further adjustment to the exchange couplings determined in Ref. [T6] . 

2. Anisotropic Exchange Model 

In this section, we provide a description of a candidate theory for the magnetism of 
Yb 2 Ti 2 7 introduced in Ref. [16]. The magnetic Yb 3+ ion has electronic configuration 
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2 Fj/ 2 , such that J = 7/2 and Lande factor gj = 8/7. The nearest neighbour distance 
between Yb 3+ ions is r nn = (a/4)\/2, where a = 10.026 A is the size of the conventional 
cubic unit cell [H]. The single ion crystal field interaction H c { is the largest magnetic 
energy scale. The form of H c { is fixed by the symmetry of the Yb 3+ environment. To 
specify H c { completely, we require the values of the six crystal field (CF) parameters. We 
use two independent sets of CF parameters determined in Refs. [H [17]. The magnetic 
Hamiltonian is H = H c { + H int , which includes both the CF and spin-spin interactions 
of the form H hlt = H cx + H dip which is composed of exchange, H cx , and long range 
magnetostatic dipolar interactions H dip = Ei> i; a,6 ^§£r I 3 ? ' J j ~ 3 ( J ? " % & )( J i ' 

with coupling D = ^ffi ~ 0.01848 K. H cx contains all nearest neighbour exchange 
interactions, Je, that respect lattice symmetries. There are four such nearest neighbour 
interactions [TH]: i?i S i ng = — Ji S i n g j> a b ( 3 i ' { 3 j ' , which couples the local 
[111] z components of J (shown in Fig. QJa)), H iso = - J iso Y^<i,j>-a,b J i ' J i> tne 
standard isotropic exchange, H pd = - J pd T,<ij>; a ,b[ J i ' 3 ) ~ 3 ( J ? ' ^ b )( J i " a 
pseudo-dipolar interaction of exchange origin and not part of H dip and, finally, Hdm = 
— cTdm X^<i j> a b ^dm ' {^i x Jj) j ^ e Dzyaloshinskii-Moriya (DM) interaction [19]. In 
all of these terms, J" denotes the angular momentum of the Yb 3+ located at lattice R" 
(FCC lattice site i and tetrahedral sublattice site a) and R"- is a unit vector directed 
along R b j - R^. 

In our recent work [16], we were able to determine the values of J e = 
{Jiso, iTising, Jpd, Jdm}- We accomplished this by fitting the diffuse paramagnetic 
neutron scattering measurements to neutron scattering patterns generated using the 
random phase approximation [16J. The fitting was done by treating the four exchange 
couplings, J7e, as free parameters, and determining their optimal values using simulated 
annealing to minimize the difference between the model and experimental scattering 
patterns, along with a term to fix 9cw to the value reported in Ref. [IT] of 6*cw = 0.75 
K [20]. The values of J e we extracted are [IB] : 

Ji sing = 0.81 K J iso = 0.22 K Jp d = -0.29 K J DM = -0.27 K (1) 

when using the CF parameters of |17j . and 

Jising = 0.76 K J iso = 0.18 K Jp d = -0.26 K J DM = -0.25 K (2) 

when using the CF parameters of [8]. The values of Jdm in both of these sets of 
interactions are large, even compared to values found in other insulating pyrochlore 
materials based on magnetic d-shell transition metal ions that display significant DM 
interactions [2U|22]. However, it has recently been shown that large DM interactions are 
allowed in Yb2Ti 2 07 due to antisymmetric contributions to the superexchange between 
Yb 3+ ions mediated via the oxygen ions at the 01 site (i.e. in the middle of each 
tetrahedron) [TTj . 
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3. Local Susceptibility 

We briefly outline the method of Gukasov and Brown [23] for determining the 
components of the site susceptibility tensor from polarized neutron scattering. In the 
presence of a magnetic field, the induced moment on sublattice a at fee site R« is given 
by M (Rj) where we assume that there is no contamination from k ^ magnetic 
structures. The magnetic Bragg intensity /(q) depends on M(q), the magnetic moment 
in reciprocal space, and is given by 



where T(q) = |aF(q)[M(q) — q(M(q) ■ q)] and both b and a are constants. Since the 
neutron intensity depends on the neutron polarization vector Pq, measurements of the 
spin-non-flip and spin-flip intensity for different Bragg peaks provide information about 
the moments on each sublattice. 

The site susceptibility Xa is defined as M°(Rj) = y^fh^ where h 13 is the applied 
magnetic field and a, (3 are components in the crystallographic frame, with the bar 
over Xa indicating the choice of the crystallographic frame. In the refinement of the 
moments from the neutron data, Xa on sublattice 1 (at the fee sites) was taken to have 
two components with all diagonal elements equal and all off-diagonal elements equal in 
the crystallographic frame. In the experiment [8] on Yb2Ti 2 07, the applied magnetic 
field was 1 T in the [110] direction, so the experiment was sensitive to components 
xf with a = 1,2. Figure [T](a) shows the tetrahedral basis with the choice of sublattice 
labelling as well as the local Ising (111) directions and the [110] orientation of the applied 
magnetic field. 

We computed the moments and hence the site susceptibility within a local mean 
field theory by computing the single ion spectra self-consistently in the presence of a 
magnetic field. Each bilinear spin operator J°J^, with spin components a and /3, in the 
interaction Hamiltonian ifi nt was written as 



and the last term on the right-hand side - the fluctuation term - was dropped to obtain a 
decoupled mean field Hamiltonian i?i n t,MF = J2i -^mt,MF(i) for each site i. The single ion 
wavefunctions \u), at site i, were obtained from Hmf = H c { + Hz + i?int,MF 0, EG] where 
Hz is the Zeeman interaction term, Hz = —^b9j J a • B by solving for the moments 
(Jf). The infinite lattice sum of the dipolar interaction was computed using Ewald 
summation [23]. From a set of randomly chosen initial moments, Hmf was diagonalized 
on each site of a cubic unit cell with periodic boundary conditions and the resulting 
spectrum, \v,i), for sites i = 1, . . . , 16 and v — 1, . . . , 2 J + 1 used to obtain a new set 



of moments (J™) = Z 1 Tr {J" exp(— /3iJ MF (i))} where Z = Tr {exp(— (3H MF (i))}. The 



process of successive diagonalization and computation of the moments was iterated until 
convergence was reached. The process was repeated for each temperature. The use of a 
single cubic unit cell for the diagonalization is more than sufficient since it was found in 




J? J? -> Jf(J?) + (J?) J? - (JfXJf) + (J? - (J?)) (J? - (J?» 
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Ref . [H] that the interaction Hamiltonian H int with parameters from either Eq. [T] or [2] 
leads to q = ordering, signifying an identical moment configuration on each tetrahedral 
primitive unit cell. 

From the moments obtained within MF theory, we computed the single ion 
susceptibility Xa for each sublattice a = 1,2,3,4 - each independent (primitive) 
tetrahedron being identical. In common with Refs. [H |9], we present components of 
the susceptibility parallel (%|i) and perpendicular (x±) to the local Ising directions. 
These components were computed by introducing rotation matrices u^f from the global 
to the local frame with along the (111) directions. Transforming \ a to the local 

frame x^f = u T u a. 5 X 1 a '> we extract Xa,\\ = X z a z ■ O ur choice of local frame is such that the 
[110] field has vanishing local y component for sublattices 1 and 2 (along the a chains) 
and vanishing local x and z components for sublattices 3 and 4 (along the /3 chains). 
Hence, for sublattices 1 and 2, with an applied [110] field, we take Xa,± = X xx and for 
sublattices 3 and 4, Xa,i. = X vv ■ 

4. Results 

Fig. DJb) shows the single ion Xa = {xbX±}> computed only from the Zeeman term 
and CF Hamiltonian using the CF parameters of Refs. [BJ [17], and the experimentally 
determined values of Xa [S]- Both sets of CF parameters give rise to a x± that fits 
the experimental data well over the whole temperature range even in the absence of 
interactions except for the lowest temperature point, T = 2 K, where the experimental 
Xj_ falls below the crystal field prediction. However, the x\\ calculations for the two sets 
of CF parameters are significantly different from one another and from the experiment. 
Whereas, the parameters of Ref. [8] give a x\\ that fits reasonably well between 20 K 
and 70 K, deviating from experiment at the highest and lowest measured temperatures, 
those of Ref. [TTJ do not fit well over most of the temperature range. The addition of 
anisotropic exchange of the form given in Eq. (OQ), does not improve the fit to experiment 
of x\\ computed using the CF parameters of Ref. [17]. Because of the poor fit to X||, even 
when anisotropic exchange is included, we do not use the CF parameters of Ref. [T7] 
any further in this work despite their fair success in reproducing paramagnetic neutron 
scattering measurements in Ref. [16] . Another point to be drawn from Fig. [Hb) is 
that there is a difference between the local susceptibilities for the a and (3 chains. In 
particular, x\\ is n °t probed for the 3 and 4 sublattices because the [110] field does not 
induce a moment in the local [111] Ising direction on these sublattices. Also, x±j which 
is probed on all four sublattices, is different for the a and (3 chains. As the temperature 
is raised, the cubic symmetry between different sublattices is partially restored. 

When interactions are included in the computation of Xa as described in Section [3j 
we obtain the results shown in Fig. |2j Figure |2^a) shows the local susceptibility, 
computed for two different models over the temperature range 2 K to 300 K, together 
with the experimental data of Ref. [8] . The inclusion of interactions does not significantly 
affect x± so ) f° r these results, we merely observe that the splitting between a and f3 
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Figure 1. (a) Tetrahedral basis of the pyrochlore structure showing the four anisotropy 
axes as red arrows and the applied magnetic field B orientation (orange arrow). 
Sublatticcs 1 and 2 make up the a chains and sublattices 3 and 4, the ft chains, 
(b) Local susceptibility (open symbols) as reported in Ref. |] in a 1 T field along 
[110], along with single ion local susceptibilities computed using the crystal fields of 
Ref. [8] (dashed and solid lines) and Ref. [17] (dot-dash line - the perpendicular part 
having been omitted as it almost coincides with the Ref. [8] crystal field result) in the 
same field, but without including the J e exchange interactions. 



chains is present once again because of the symmetry breaking field. The first model 
to be compared with the experimental x\\ data is one with long ranged dipoles and 
purely isotropic exchange (J7i s i ng = J7pd = Jdm = 0), with a coupling <J iso = —0.06 K 
constrained by the Curie- Weiss temperature which is taken to be #cw = 0.75 K, within 
the experimental error margins given above [16] (blue dot-dash line). The second model 
is the candidate anisotropic exchange model described in Section [2] and Ref. [16] (blue 
dashed line). The x± results (red lines) do not strongly distinguish between the two 
models. At low temperatures, T < 30 K, where the effect of interactions produces 
significant deviations of x\\ from the noninteracting case, the isotropic exchange model 
does not fit the experimental data as well as the anisotropic exchange model, except at 
the lowest temperature point (2 K). However, neither model captures the 2 K point for 



Fig. H£b) shows once again the experimental data and the anisotropic model 
predictions for Xa as well as Xa computed using the molecular field method of Ref. [8] . 
The fit in Ref. [8] of x\\ is only provided down to approximately 3.5 K in that work (and 
no results for x± for the single ion mean field case are reported in that reference). For the 
temperature range over which the models can be compared, our anisotropic exchange 
model describes the experimental data as well as the model of Ref. [8]. However, our 
model achieves this while correctly preserving lattice symmetries, unlike the model of 
Ref. [8], and therefore corresponds to a physically acceptable microscopic model. We 
discuss further the comparison of our results with those of Refs. [HI Ej in the Appendix. 

At temperatures T ^> # cw , interactions are expected to have little effect on 
Xa which is borne out by our calculations. Therefore the differences between the 



X±- 
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Figure 2. (a): \a as reported by [5] (open squares and circles). Also shown is \a, for 
all four magnetically incquivalent sublattices, computed within local mean field theory 
with our anisotropic exchange model with long range dipoles obtained from neutron 
scattering data in Ref. [16] . The temperature range is 2 K to 300 K. Also shown is 
the prediction of a model with isotropic exchange interactions and long ranged dipoles. 
(b): A comparison of our anisotropic exchange model with the prediction of the model 
of Ref. [8]. 

Xa calculations and experimental results at high temperatures (T > 100 K) can be 
attributed to an incorrect spectral decomposition of the excited crystal field states [8lfT7] . 

5. Discussion 

We have investigated the local susceptibility computed from the anisotropic exchange 
model described in Section [2] for Yb2Ti 2 07. The model was obtained from an 
independent previous fit to the paramagnetic neutron scattering [16] . The calculations 
presented here fit the Xa experimental data well with no adjustable parameters over 
much of the experimental temperature range (2 K to 250) K. There are, however, 
discrepancies in the fit above approximately 100 K, presumably due to the fact that 
the available CF parameters do not capture the correct wavefunctions for the excited 
CF energy levels. The fit is also imperfect at the lowest observed temperature of 2 K 
where the experimental Xa susceptibility flattens out while it does not in the mean field 
treatment of the Hamiltonian H = H c f + H mt + H%. 

At high temperatures T 3> #cw; Xa is insensitive to the interactions and the quality 
of fit is determined by the CF Hamiltonian and the Zeeman term. Fig. [T] illustrates that 
the crystal field provides an adequate fit to x\\ from 100 K down to about 30 K and 
to x± over t ne entire range of temperatures. Figure [2] shows that the quality of fit 
for T > 30 K does not change with the inclusion of interactions. Therefore, the only 
nontrivial test of our model comes from the five lowest temperature points (T = 2 K, 
5 K, 10 K, 15 K, 20 K) in x\\ where the noninteracting Hamiltonian, H = H c f + Hz, 
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gives a \\\ that falls significantly below the experimental result (see Fig. DJb)). The 
anisotropic exchange model determined in Ref. [16] provides a good description of 
whereas isotropic exchange interactions do not - this is the main conclusion of this work. 
Although the relevant test of our model comes from only a few (five) data points, it 
achieves the goals of providing a consistency check of our model determined previously 
on the basis of diffuse paramagnetic neutron scattering [16] while also establishing that 
anisotropy in the exchange is indeed an intrinsic ingredient of the spin Hamiltonian of 
Yb 2 Ti 2 7 . 

Refs. [8] and [9] both extract exchange parameters from the local susceptibility data, 
and present fits to the local susceptibility that result from these exchange parameters. 
As stated above, we find that this method for extracting exchange parameters is under 
constrained. In the Appendix, we describe in some detail the approaches taken in these 
two papers and compare with our own results. 

6. Conclusion 

In this work we have shown that the anisotropic exchange model extracted from neutron 
scattering measurements in a previous work [16] produces a good fit to both local and 
bulk susceptibility measurements on the rare earth pyrochlore material Yb2Ti 2 07 for 
temperatures in the range T < 100 K. This model has four nonadjustable and a priori 
prescribed bilinear exchange couplings which were extracted from neutron scattering 
data [16]. The fit of our model to the experimental measurements of \ a and the bulk 
susceptibility provides further independent evidence for the correctness of our anisotropic 
exchange model |16j . 

We also found that the origin of the discrepancy between the experiment and our 
model at high temperatures, and indirectly, a similar failure to fit the local susceptibility 
Xa also at T > 10 2 K in Ref. [8], lies in the excited CF states obtained from the CF 
parameters of Refs. [SI HZ]- Further experimental work would be required to probe more 
accurately the spectral decomposition of the excited states of the crystal field in order 
to better determine their structure. That said, we do not expect that the corrections 
to the spectral decomposition of those excited states would, as a byproduct, necessitate 
revisiting the value of the anisotropic exchange parameters J e . We are thus confident, 
combining the results presented in this paper and those of Ref. [16] , that we have in hand 
the correct microscopic Hamiltonian of Yb2Ti 2 0y that includes all symmetry-allowed 
anisotropic bilinear nearest-neighbour couplings. One may therefore aim to proceed to 
unravel the nature of the low temperature (T < 2 K) state of this exotic material. In 
this context, it is worth commenting on the validity of the mean field theory treatement 
of our model at low temperatures (T < 2 K) and/or in weak field. Mean field theory 
predicts a transition to a q = structure at approximately 1.1 K [16]. This contrasts 
with the experimental T c = 240 mK [12]. As usual, mean field calculations neglect the 
effects of both thermal and quantum fluctuations. These could be fairly large in this 
system and greatly decrease the critical temperature. The role of thermal and quantum 
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fluctuations in the proposed Hamiltonian H = H mt + H c { for Yb 2 Ti 2 07 remains to 
be investigated. To the best of our knowledge, and notwithstanding the difficulties 
with mean field theory, the sharp specific heat feature signaling a phase transition in 
Yb2Ti 2 07 has only been reported in powder samples. It would be useful to confirm the 
existence of this feature in single crystals. 

Two interesting observations that arise from this work are that (i) the isotropic 
exchange J lso is ferromagnetic (positive) in Yb2Ti 2 07 while it is found to be 
antiferromagnetic (negative) for the other ^2^207 compounds and (ii) there appears to 
exist large antisymmetric anisotropic exchange, described above in terms of the form of 
a Dzyaloshinskii-Moriya (DM) interaction. Very recent theoretical work [H] has found 
that large anisotropic exchange, including large DM interactions, can in principle arise 
in Yb 2 Ti 2 07. It would be interesting to investigate this issue further. 

The conclusion of this work, that the magnetic exchange interactions in Yb2Ti2C>7 
are anisotropic, hints at the possibility of significant multipolar interactions [TUJ, [25] in 
this material and, perhaps, in other rare earth pyrochlores. That said, since the crystal 
field gap is much larger in Yb2Ti 2 07 than the scale of the interactions, anisotropic 
bilinear couplings between the Jj operators as as well as multipolar couplings between 
the Yb 3+ moments manifest themselves, at low energies, in the form of generally 
anisotropic bilinear exchange couplings between effective pseudospins one-half [TT] . 
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Appendix A. Previous works on local susceptibility in Yb 2 Ti 2 7 

In this Appendix we discuss the anisotropic exchange fits of Refs. [Sj, for Yb2Ti 2 07, and 
[9] for various rare earth pyrochlore titanates. 

The fit to Xa presented in Ref. [8] and reproduced in Fig. [2b is based on a self- 
consistent mean field theory that considers a single sublattice. The moment (J), at a 
temperature T, was computed from a mean field Hamiltonian consisting of the crystal 
field, a Zeeman term, and a mean field term arising from the interactions of the form 
—gj^B^ai^a) j where a specifies the components || and _L, the directions parallel and 
perpendicular to the local [111] direction. Solving for the moment, the two couplings 
A || and Aj_ were determined from a fit to Xa- Since this model fails to take into account 
the sublattice structure of the pyrochlore lattice, the ensuing fit is unlikely to correctly 
reflect the underlying microscopic couplings in Yb2Ti207. 

The model of Ref. [9] includes the crystal field anisotropy as well as nearest 
neighbour exchange couplings and the long range dipole interaction (treated via an 
Ewald summation). The model correctly handles the sublattice structure of the 
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pyrochlore lattice. The model includes, at the outset, three exchange couplings which 

\™ 2 and A J 



respect the lattice symmetries and which are denoted: A^ 1; \¥ 2 an d • These can be 



expressed in terms of the couplings J c as follows: 

A_l x = jTiTlsing "I" c7iso + cTpd C^-l) 
A^ 2 = <^iso + |<^pd (A-2) 
A||^ = — fiTlsing + >J'mo ~~ 2j7pd (A-3) 

This set of couplings does not include the Dzyaloshinskii-Moriya coupling, J7dm- 
Furthermore, in the fits to Xa, Ref- [H] imposes the constraint = \ K f 2 - While 
this choice most likely circumvents the problem of having an underconstrained set of 
couplings to fit the Xa data, the reduction of the parameter space to a particular two 
dimensional surface within a four dimensional space of couplings appears to us to be 
somewhat arbitrary. Indeed, the couplings we have obtained from neutron scattering 
data and which we used in our fit to Xa do not lie within the parameter space considered 
in Ref. [9]. For example, Jbm — —0.25 K for our model with crystal field parameters 
taken from Ref. [8] (see Eq. [2]), which is of the same order of magnitude as the other 
three couplings. Reference [9] makes use of their model to study the exchange couplings 
of a number of rare earth pyrochlore titanates finding a significant degree of anisotropy 
in the exchange in a number of cases. For the case of Yb 2 Ti 2 07, which is relevant to this 
article, they mention couplings that are much larger than those of the other materials 
considered but do not report these couplings. 
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